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We present the surface worm algorithm (SWA) which is a generalization 
of the Prokof 'ev Svistunov worm algorithm to perform the simulation of 
the dual representation (surfaces and loops) of Abelian gauge-Higgs models 
on a lattice. We compare the SWA to a local Metropolis update in the dual 
representation and show that the SWA outperforms the local update for a 
wide range of parameters. 
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1. Introduction 



The complex fermion determinant at finite chemical potential has slowed 
down the progress in the exploration of the QCD phase diagram using Lat- 
tice QCD. Among the different techniques to deal with the sign problem (see 
e.g. [I]), the dual representation is a powerful method which can help us to 
solve the sign problem without making any approximation of the partition 

t" — sum as in other methods. Before approaching the ultimate goal of finding a 

dual representation for non- Abelian gauge theories with only positive prob- 
ability weights, which is a rather involved task and has not been achieved 
yet, it is advisable to explore and understand the method in simpler models, 
such as Abelian theories coupled to scalar fields [2J |3] that we study here. 

y—{ Once the partition sum with real and positive probability weights is 

found, the next step is to choose the most efficient algorithm to save com- 
puter time. In the case of only matter fields or spins, the worm algorithm 
[3] constitutes one of the most suitable and efficient methods [5] to deal with 
the constrained degrees of freedom of the system, i.e. with loops. In this 
article we present an extension of the worm algorithm (SWA) [6] to perform 
the simulation of the U(l) gauge-Higgs model where loops and surfaces are 
the dual variables. We assess the performance of the SWA in comparison to 
a local Metropolis update (LMA). The analysis of the physics of this model 
will be presented elsewhere [7j. 
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2. Dual representation of an Abelian gauge-Higgs model 

The action of the U(l) gauge-Higgs model on the lattice is given by 

X V<p 

Sh = E [k\4> x \ 2 + M<Px\ 4 } - E [4>* x u x ,u4>x+i> + 4>* x u*_c iV (j) x+ p] (i) 

where U X)V = e lAv G U(l), A u G [— vr,7r] are the link variables, 5g is the 
usual plaquette action, with U x>up = U XiU U x+ v tP U* + ~ u U* p the plaquette 
variable. The Higgs fields 4> x G C in the Higgs action Sh live on the sites 
of the lattice, k is a mass parameter and A is the quartic coupling. 

Here we outline the general strategy for the derivation of the dual represen- 
tation (for the details see the appendix in [6]). The general steps are: Write 
the Boltzmann weight in a factorized form and expand the exponentials for 
individual plaquettes and links. A single nearest neighbor term turns into 

A single plaquette term leads to 

pu X! vU x+i> ,pU* +fi u* _ sr^ P x ' vp i>r j j j j* r j* iPx,u P 

e ~~ 2_^ r, ! I x,i>U x+v,p<~> x +p,v u x,p\ 

Px,up- 

Px,up 

After integrating out the U(l) variables, the new form of the partition sum 
depends only on the dual variables: The constrained link occupation number 
lx,v G (—00, +00), the unconstrained link occupation number l X)V G [0, +00) 
and the constrained plaquette occupation number p XjUp G (—00, +00). The 
new form of the partition sum is 

Z ex Y, yV\p,l,l]C s [l]C L [p,l], (2) 

{P,U} 

where the new degrees of freedom are the dual variables /, I and p and 
Yltp 1 1} denotes the sum over all their configurations. W[p, 1, 1} is a positive 
weight factor. Furthermore, constraints appear that force the total sum of 
the occupation numbers to vanish at every site and link: C$[l] is the site 
constraint which forces the total matter flux to vanish at every site and 
gives rise to loops, 
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The link constraint Cl\p, I] gives rise to gauge surfaces, 

\/x, v : ^2 \p x>up - p x -p,up] ~ ^2 ^P X 'P V ~ Px-p'P"} + 1*> V = ° • 

\p:u<p p:v>p J 

3. Monte Carlo simulation 

To perform the Monte Carlo simulation of the system we developed the SWA 
and we compared its performance against a local update (LMA) [3j [6] . The 
LMA consists of: 

• A sweep of the unconstrained variables I rising or lowering their oc- 
cupation number by one unit. 

• "Plaquette update" : It consists of increasing or decreasing a plaquette 
occupation number p x ,vp and the link fluxes l X;(T at the edges of p x ,up 
by ±1 as illustrated in Fig. [T] The change of p x ,up by ±1 is indicated 
by the signs + or — , while for the flux variables we use a dashed line 
to indicate a decrease by —1 and a full line for an increase by +1. 

• "Cube update" : The plaquettes of 3-cubes of our 4d lattice are changed 
according to one of the two patterns illustrated in Fig. [2] Although 
the plaquette update is enough to satisfy ergodicity, the cube update 
helps for decorrelation in the region of parameters where the link I ac- 
ceptance rate is low and the system is dominated by closed surfaces. 

A full sweep consists of visiting the 4V4 links, 6V4 plaquettes and 4V4 3- 
cubes, offering one of the changes mentioned above and accepting them 
with the Metropolis probability computed from the local weight factors. 



Fig. 1. Plaquette update: A plaquette occupation number is changed by +1 (lhs. 
plot) or —1 (rhs.) and the links of the plaquette are changed simultaneously. The 
full line indicates an increase by +1 and a dashed line a decrease by — 1. The 
directions 1 < v\ < V2 < 4 indicate the plane of the plaquette. 

Instead of the plaquette and cube updates we can use the worm algorithm. 
The SWA (see [6] for a detailed description) is constructed by breaking up 
the plaquette update into smaller building blocks called "segments" (ex- 
amples are shown in Fig. K3j) used to grow surfaces on which the flux and 
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Fig. 2. Cube update: Here we show the changes in the plaquette occupation num- 
bers. The edges of the 3-cube are parallel to the directions 1 < V\ < v 2 < v?, < 4. 

plaquette variables are changed. In the SWA the constraints are temporar- 
ily violated at a link Ly, the head of the worm, and the two sites at its 
endpoints. The admissible configurations are produced using 3 steps: The 
worms starts by changing the flux by ±1 at a randomly chosen link Lq (step 
1 in Fig. HI). Lq becomes the head of the worm Ly. The defect at Ly is then 
propagated through the lattice by attaching segments, which are chosen in 
such a way that the constraints are always obeyed (2 in Fig. El). The defect 
is propagated through the lattice until the worm decides to end with the 
insertion of another unit of link flux at Ly (3 in Fig. HI). A full sweep with 
the SWA consists of V4 worms. 




Fig. 3. Examples of positive (lhs.) and negative segments (rhs.) in the ^-z/2-plane 
(v\ < v-i). The plaquette occupation numbers are changed as indicated by the 
signs. The full (dashed) links are changed by +1 (—1). The empty link shows 
where the segment is attached to the worm and the dotted link is the new position 
of the link Ly where the constraints are violated. 



<A 



Fig. 4. Illustration of the worm algorithm. See text for an explanation. 



4. Numerical analysis 

For the comparison of both algorithms we analyzed the bulk observables 
(and their fluctuations): Up which is the derivative wrt. f3 and \cp\ 2 (deriva- 
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Fig. 5. Observables as a function of /? for different parameters and volumes. We 
compare results from three algorithms: The conventional approach (asterisks), the 
SWA (circles) and the LMA (crosses). 

tive wrt. k). First we checked the correctness of the SWA comparing the 
results for different lattices sizes and parameters. For example, the upper 
plot of Fig. [5] shows Up as a function of (3 for k = 4, A = 1.5 on a lattice 
of size 4 4 . The lower plot shows (|<^| 2 ) for k = 8 and A = 1 on a 10 lat- 
tice. In both cases we used 10 6 equilibration sweeps, 10 6 measurements and 
10 sweeps for decorrelation between measurements. We observe very good 
agreement among the different algorithms. 

In order to obtain a measure of computational effort, we compared the nor- 
malized autocorrelation time r as defined in [6] of the SWA and LMA for 
different volumes and parameters. We concluded that, the SWA outper- 
forms the local update near a phase transition and if the acceptance rate 
of the link variable I is not very low (eg. lhs. of Fig. k3|). On the other 
hand, when the links become expensive the worm algorithm has difficulties 
to efficiently sample the system (as can be observed on the rhs. of Fig. pi r 
for Up is larger for the SWA than for the LMA). But this can be overcome 
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Fig. 6. Normalized autocorrelation times r for 2 different set of parameters. Left: 
parameters close to a first order phase transition. Right: low acceptance rate of 
the variable I. Both simulations correspond to a 16 4 lattice. Data taken from [5J. 

by offering a sweep of cube updates or a worm made of only plaquettes as 
described in |2]. 
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